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We complete the introduction of Yano-Koonin-Podgoretskii (YKP) parametriza- 
tion of the correlation function by deriving the corrections to the correlation radii 
which arise from releasing the popular smoothness approximation and approxima- 
tion of setting the pair momentum on-shell. We investigate the definition range of 
this parametrization and find kinematic regions in which the YKP parametriza- 
tion is inapplicable. These problems disappear if the newly proposed Modified 
Yano-Koonin-Podgoretskii parametrization is used. We then focus on the physical 
interpretation of the correlation radii obtained in the different parametrizations. 
While the extraction of the longitudinal source expansion from the YK rapidity 
is found to be rather robust against variations of the source density profiles, the 
extracted emission duration is quite sensitive to such variations and cannot be 
reliably extracted except for the case of extremely long source lifetimes. 
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1 Introduction 

In recent years Bose-Einstein interferometry [1, 2] for ultrarelativistic heavy ion colli- 
sions has been developed into a powerful analysis tool [3, 4]. It turned out that not only 
the sizes of boson emitting sources are measurable via the magnitude of the so-called 
correlation radii, which are the parameters of a Gaussian fit to the observed correla- 
tion function, but also the source dynamics leaves its fingerprint in their dependence 
on the average momentum of the pair [5, 6, 7, 8, 9]. On the other hand it became also 
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clear that the temporal and spatial size parameters of the emitter are not all separately 
accessible; only certain combinations of the space-time variances of the source are mea- 
surable [10, 11, 12]. This knowledge is crucial for an understanding of the usefulness 
of various Gaussian paramctrizations of the correlator. By an appropriate choice of 
the parametrization one can ensure that the resulting width parameters reflect specific 
combinations of the space-time variances which allow for a particularly straightforward 
physical interpretation in terms of specific source features. 

It has recently been argued that for sources with dominant longitudinal expansion 
the Yano-Koonin-Podgoretskii (YKP) parametrization of the correlator [13, 14, 15] pro- 
vides a more direct interpretation of the measurement than the "traditional" Cartesian 
one [14, 10, 11]. It was stated that the YKP parametrization allows for direct ob- 
servation of the longitudinal expansion of the source and provides a more accurate 
measurement of the emission duration. These arguments were supported by a study 
involving a model which was believed to reproduce the data and to include all relevant 
features admitted by experimental observations [12, 16]. Of course, strictly speaking 
the "amount" of information gained from a measurement is independent of the chosen 
parametrization of the correlator (as long as it is complete); the extraction of certain 
quantities of interest appears, however, easier when using the YKP one. The mathemat- 
ical equivalence of the different parametrizations leads, on the other hand, to relations 
between their respective correlation radii which must be satisfied if the analysis was 
done correctly. The corresponding relations between the Cartesian and YKP radius 
parameters were given in [12, 16]. 

In this paper we complete the program of introducing the YKP parametrization. 
Corrections to the correlation radii due to the off-shcllness of the pair momentum K 
and due to the so-called smoothness approximation are given up to 0{q 2 ). (The corre- 
sponding corrections for the Cartesian parametrization were derived in [17].) Moreover, 
we study the definition range of this parametrization. We find that kinematic regions 
can exist where it is not defined. A small modification in the formulation of YKP 
parametrization is shown to avoid these problems. The cost for this remedy is a slightly 
less straightforward interpretation of the resulting radius parameters. 

With these improved results at hand we can once more study the question of the 
interpretation of the measured correlation radii in the different parametrizations. We 
focus in particular on the extraction of the source lifetime (i.e. the emission duration). 
We find that that the original optimistic statements in this context [12, 16] require 
qualification - the lifetime measurement turns out to be considerably more model- 
dependent than originally anticipated. We illustrate this with some simple examples in 
Section 3. 

More technical items are deferred into the appendices. In Appendix A the definition 
range of YKP parametrization is studied. Relations between the correlation radii of 
different paramctrizations used in this paper are given in Appendix B. Finally, the 
mentioned 0(q 2 ) corrections to the YKP radius parameters are derived in Appendix C. 
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2 Measuring the source 

Given the emission function S(x, K) (Wigner phase-space density of the source) and 
assuming completely chaotic (independent) production of the particles [23], the corre- 
lation function can be calculated via [18, 19, 20, 21, 22] 

C (gK) = l I \Jd^e^S( X ,K)f 

Here K and q are the average and relative 4-momentum of the particle pair: 

K = ^(pi+P2), q=Pi-p2- (2) 

Equation (1) requires the knowledge of the full quantum mechanical structure of the 
source: the momentum K in the numerator is off-shell. In practice, however, when 
evaluating the correlations it is put to its on-shell value by setting 

K°^E K = VK 2 + m 2 . (3) 

The formalism also simplifies considerably if both momenta in the denominator of (1) are 
identified with K (the so-called smoothness approximation) . The applied simplifications 
are expected not to be too serious because the interesting region of correlations lies at 
small values of q where the above approximations are good. Nevertheless, corrections 
due to the above simplifications can be (and will be) derived, which support these 
superficial statements in a more rigorous way. 

To summarize, the above two approximations lead to the relation 

with K taken on-shell according to (3). 

It is convenient to parameterize the correlation function in Gaussian form: 

C(q,K)-l = e X p(-Q^B^(K)). (5) 

Expressions for the elements of the (symmetric) matrix B are obtained by expanding 
both (4) and (5) up to second order in q and identifying the coefficients of both series. 
This leads to 

B^u = {x^x v ) , (6) 

where 

Xn x^ (Xfi) j (7) 
and (. . .) denotes the average weighted by the emission function: 



(8) 
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One must, however, keep in mind that correlations of two physical particles are 
measured and therefore q and K fulfill the on-shell constraint 

q-K = 0, (9) 

which reduces the number of independent g-components to three. Due to this only 
certain combinations of the B^s are experimentally accessible. 

In what follows we focus entirely on a class of azimuthally symmetric sources repre- 
senting the fireballs generated in central collisions. It is common to orient the coordinate 
system in such a way that the z-axis (longitudinal, I) lies in beam direction and the 
x-axis (outward, o) in the direction of the transverse component of the pair momentum. 
The remaining y-direction is then denoted as sideward (s). For azimuthally symmetric 
sources there exists then a symmetry under the transformation y — > — y, leading to 
only seven non-zero components of B of which only four combinations are measurable 
[15]. (Note that no such symmetry is present in the outward direction as long as the 
directions parallel and antiparallcl to K± are not equivalent.) 

In this coordinate frame q s decouples from the equation (9). When employing 
this relation in different ways one arrives at different Gaussian parametrizations of 
the correlator. 

2.1 Cartesian parametrization 

Resolving the on-shell constraint (9) as 



qt = fi±q + fiiqi , (10) 

where /3j = K l / Ek and qt is the zeroth component of relative momentum, and inserting 
(10) in (5) leads to the Cartesian parametrization of the correlator [14, 10, 11] 

C(q, K) - 1 = C M-q 2 s R 2 s (K) - q 2 Rl(K) - q 2 R 2 (K) - 2q q l R 2 ol {K)] . (11) 

For the correlation radii one obtains 

Rl = ((x-/? ± t) 2 }, (12a) 

R 2 S = (y 2 >, (12b) 

R 2 = ((z -Pit) 2 ), (12c) 

R% = ((x-/3 ± t)(z-At)>. (12d) 



Due to the on-shell constraint (10) spatial and temporal source sizes are mixed in 
all observable radius parameters except for R 2 S . This makes the interpretation of a 
measurement rather involved since the spatial and temporal contributions cannot be 
resolved in a model-independent manner. The situation simplifies for the longitudinal 
radius R 2 if the frame in which analysis is performed is boosted longitudinally such 
that K t = A = (Longitudinally Co-Moving System, LCMS) [24]. Then 



Rf - <2 2 >, 

Rll = ((x-#lt)z) 



(13c) 
(13d) 
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So the longitudinal dimension, in the longitudinal rest frame of the pion pair, is mea- 
sured directly. The cross-term i? 2 ; may obviously take both positive and negative values. 
The superscript should not confuse - it just expresses the bilinear character (and the 
dimensionality) of this parameter. 

Access to the lifetime is provided by the difference [24] 

R 2 diS = Rl-R 2 s =Pl (t 2 ) - 2/? ± (it) + (x 2 - f) . (14) 

Of course, one has to assume that i? 2 iff is dominated by the first term on the r.h.s., at 
least at higher K± where the suppression by (3± is small. In Section 3 we investigate 
this assumption in more detail. 

Note finally that corrections to the correlation radii given by (12) due to the ap- 
proximation (4) were derived in [17] - they are listed in Appendix C. 

2.2 Yano-Koonin-Podgoretskii parametrization 

When one wants to measure the source using Bose-Einstein interferometry the problem 
of the "right" frame for the measurement arises. The appropriate choice of frame may 
even depend on the pair momentum. Only a part of the whole (expanding) fireball 
size is reflected in the correlation radii, namely that part contributing to production of 
particles of a given momentum [5] . In the following this particular part of the source 
will be called effective source. 

Since the measurements are done in order to obtain information on the local phase- 
space characteristics (density, energy density, temperature etc.) of the emitter one 
should expect that the appropriate frame for the measurement of the (effective) source 
sizes is its own rest frame. Unfortunately, this frame cannot be directly determined in 
the experiment. However, in [14] Podgoretskii showed for non-expanding sources that 
the cross-term of the correlator parametrization is a direct consequence of measuring out 
of the source rest frame. The issue becomes a bit more involved when expanding sources 
are treated, because then expansion is present also in the effective source. However, 
since the velocity difference between its edges is not very big, one can still (at least 
intuitively) introduce something like an average or effective source velocity and try to 
express the cross-term with its help. 

This can be to a large part achieved by resolving the on-shcll constraint (9) as 

1 # (in 

q = -fi-qt - -z-qi ■ (15) 

and rewriting (5) in the so-called Yano-Koonin-Podgoretskii (YKP) form [15, 12]: 

C(q,K) - 1 = eM-qlRl(K) - (qf-q 2 t )R 2 (K) - (R 2 (K)+R 2 (K))(q ■ U(K)) 2 }. (16) 

Here the three correlation radii are manifestly invariant under longitudinal boosts, and 
the cross-term has been absorbed into a term involving the Yano-Koonin (YK) velocity 
v YK which is the longitudinal component of the 4- velocity U: 



[7 = 7(1,0,0,0, 7 = pr— — 2 — • d~) 



G 



Boris Tomasik, Ulrich Heinz 



Note also that in the Yano-Koonin-Podgoretskii paramctrization the components of the 



transverse momentum difference occur only in the combination 

Q± = VW+tf- (18) 

For the YKP correlation radii and the YK velocity one can infer expressions similar 
to (12) [12, 16]: 

Rl = A-v YK C, (19a) 

R\ = (y 2 )=i? 2 , (19b) 

R\ = B-v YK C, (19c) 

- - ^(>-H^) 2 )- 

Here we introduced the shorthands 

A = (t 2 )- A(xt) + i^(x 2 -y 2 ), (20a) 

B = (z 2 )-2A(xz) + |-(x 2 -y 2 ), (20b) 

C = (tz)--l(xz)- A(xt) + A(x 2 -y 2 ). (20c) 



Note that these shorthands would appear as parameters of a (pseudo-) Cartesian para- 
metrization with the same g-components as used in YKP [9] . 

Although the YKP radii are invariant under longitudinal boosts, their interpretation 
is connected with the frame where v YK = 0, the so-called Yano-Koonin frame. In 
this frame C vanishes. From (19a) and (20a) one concludes that Rq (and only this 
parameter) is sensitive to the lifetime (t 2 ) of the source in this frame. This has led to 
some optimism in Refs. [12, 16] because measuring the lifetime in this way seems to be 
much easier than via i? 2 iff (cf. eq. (14)), due to the missing factor 0]_ which may be 
small. (Note that in the YK frame R\ is = (3±A.) However, as already mentioned in 
the discussion below (14), both extraction procedures strongly rely on the assumption 
that the last two terms in (20a) are small. For small K± this is not so clear, due to the 
small factors in the denominator: for example, although at K± = we generally 
expect (x 2 — y 2 ) to vanish due to azimuthal symmetry, if divided by {5\ the limit for 
K± — ► can be a not-zero quantity [16]. The true behaviour of A for small but nonzero 
K± can thus only be clarified by model studies. For large K±, on the other hand, 
different models give different predictions for (x 2 — y 2 ) (see next Section). Unless (t 2 ) 
becomes very large (much larger than the geometric contributions (x 2 ), (y 2 )), which 
may happen during compound nucleus formation at low energies [26] but is unlikely 
for the rapidly disintegrating sources formed in relativistic heavy ion collisions, the 
extraction of the lifetime thus appears to be quite model-dependent. This reduces the 
power of Rq (and R^ iff as well) for measuring the emission duration. In fact, Rq may 
even become negative due to the discussed "correction term" [3] . 
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On first sight, since (20b) has a similar structure as (20a), one might conclude that 
the same problems must appear when extracting the longitudinal source size from R 2 . 
Fortunately, here the situation is better since now the sub-leading terms are multiplied 
with Pi and f3f, respectively. In the YK frame this velocity is usually small (see below). 
Although this is again a model-dependent statement, it appears to be valid for all 
models which are in agreement with the data. Thus R 2 is a good estimator for the 
longitudinal dimension of the source measured in YK frame. 



2.3 Modified Yano-Koonin-Podgoretskii parametrization 

So far, nothing has been said about the definition range and the resulting applicability 
of YKP parametrization. Clearly, it is not defined for K± = (3± = as the on- 
shell constraint (15) is singular at this point. However, a more serious problem is 
that the argument of the square root in (19d) can become negative (see Appendix A). 
Then the YKP parametrization is ill- defined (the YKP parameters become complex 
and thus unphysical). We hasten to stress that the applicability region of the YKP 
parametrization cannot be checked directly by the experiment. It is therefore strongly 
recommended to cross-check the correlation radii obtained from a YKP fit via the 
relations (50)/(51) and/or (48)/(49). 

As shown in Appendix A, the failure of the YKP parametrization is due to the use 
of q± instead of q s as an independent variable. In other words, the elimination of q via 
(15) is not done completely. This leads to the appearance of contributions from (f 2 ) / 
in (20). One can, however, eliminate q D completely without giving up the longitudinal 
boost-invariance of the fitted radius parameters by writing the correlator in the form 

C(q,K)-l = cxp[-q 2 R' 2 (K)-(qf-q 2 )R\f(K) 

-(R' 2 (K) + R' 2 (K))(q-U'(K)) 2 }, (21) 

with 

U' = 7 '(1, 0, 0, <J , V = 1 . (22) 

This will be called the Modified Yano-Koonin-Podgoretskii (YKP') parametrization. For 
its correlation radii one finds 

R'o = A'-v'^C, (23a) 
R'l = (f)=Ri=R 2 , (23b) 
= B '-v'C', (23c) 



where 



A' = <t 2 )-A(xt) + ^(x 2 )=A+^<y 2 >, (24a) 
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(24b) 



(24c) 



The interpretation of R' 2 and i?y 2 is again easiest in the frame where v' YK = 0, the 
"Modified Yano-Koonin (YK') frame" . It is not as straightforward as in case of the 
original YKP radii. This is the price for a parametrization which is defined everywhere. 
While in YKP parametrization one could hope that the last "correction terms" of 
(20a), (20b), and (20c) are small, now only the (large) term (x 2 ) appears instead of the 
(smaller) difference (x 2 — y 2 ). Therefore, the deviations of R' Q 2 and R' 2 from (t 2 ) and 
(z 2 ), respectively, are larger than those of Rq and R 2 . Still, for the measurement of 
the longitudinal size (z 2 ) the situation may not be bad because (x 2 ) in (24b) is again 
multiplied by (3?, which is small in the YK' frame, exactly like discussed above in case 
of A measurement of (t 2 ) via R' 2 , on the other hand, is definitely impossible. 

A good feature which is preserved in the Modified YKP parametrization is that the 
velocity v' YK still maps to a very good approximation the velocity of the effective source. 
This will be demonstrated in the next Section. 

Finally, we note that while the relative momentum components used in the original 
YKP parametrization satisfy the inequality 



which means that the data points never fill the whole three-dimensional g-space, no 
such restriction exists for q s which is used in the modified YKP parametrization. This 
should help to avoid certain technical problems in the fitting procedure which can occur 
with the YKP parametrization. 

Cross-relations Since all parametrizations discussed so far are just different repre- 
sentations of the same correlation function, there must exist relations between their 
correlation radii. They are listed in Appendix B. 

These relations provide a powerful check of the correctness of results obtained from 
fitting the measured correlation function with Gaussian parametrizations; if Cartesian 
and YKP (and eventually Modified YKP) fits are performed independently, the resulting 
correlation radii must fulfill the corresponding cross-check relations. Such control is 
strongly recommended in order to avoid the pitfalls related to the possible non-existence 
of the YKP parametrization for the data under study (which may not show up clearly 
in the fitting process but might cause it to converge to a wrong result). 

Off-shell corrections Before closing this section we study the corrections to the 
YKP and Modified YKP correlation radii connected with the approximations made 
in (4). They are easily inferred by combining the known expressions (55) and (56) 
according to (51) and (53). The resulting corrections are given in Appendix C. In 
the case of the Modified YKP parametrization an interesting phenomenon is observed: 



1 Pi 
- q ° = 0l qt -pl 



(25) 
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while the correction terms accounting for the release of the smoothness approximation 
Pi(pi)Pi(P2) ~ [-Pi(^)] 2 m the denominator (where Pi(p) = J d 4 xS(x,p) is the in- 
variant single-particle spectrum) and those accounting for the release of the on-shell 
approximation K° w Ek in the numerator diverge separately in the limit K± — > 0, this 
singularity cancels in their sum. For the YKP paramctrization the singularities already 
cancel at the level of the individual corrections. The reason for this difference is the use 
of q± in the YKP parametrization (instead of q s as in the modified one), leading to the 
appearance of i? 2 iff in (51) (instead of i? 2 in (53)). 



3 Examples 

In this Section we illustrate some statements made in the previous Section with the 
help of a model for a longitudinally and transversely expanding fireball which is locally 
thermalized. This is neither an extensive model study nor an analysis of real data. 
Therefore we will stay on a rather superficial level and won't go very deeply into the 
model structure. 

Our model (in fact it is a class of models) is expressed by the emission function 
[8, 9, 16] 

S(x, K) d 4 x = —^M ± cosh(r - 77) exp ' K ' U ^ 
(27r) d 

xGWex P f-fc# 



exp ( — ^ T ) dnrdrdip. (26) 



V2tt At 2 V 2 A t 



The momentum is parameterized here in terms of transverse mass M± = y 1 K\ + to 2 , 
rapidity Y, and the transverse momentum K± such that 

K = (M± coshF, K ± , 0, Mi sinhF) . (27) 

As space-time coordinates in the transverse plane we use the usual polar coordinates 
r and <p, and the remaining two directions are parameterized by longitudinal proper 
time t = Vt 2 — z 2 and space-time rapidity r/ = 0.51n[(t + z)/(t — z)]. The ^-profile is 
chosen to be Gaussian with width A77 and the center at 770. The freeze-out times are 
distributed by a Gaussian of width At centered around the average freeze-out time to- 
We will study models with two different transverse profiles: a Gaussian one 

and a box-shaped one 

G(r) = 6{R B - r) . (29) 
The dynamics of the source is implemented via the velocity field u(x), 

u(x) = (coshrjt cosh?7, cos 92 sinh?7i, sirn^ sinh?7 t , cosh?7 t sinhr/) , (30) 
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Table 1 Values of model parameters used in the calculation 

box Gauss 

temperature T 0.1 GeV 

transverse flow scaling parameter 77/ 0.57 0.64 

geometric transverse radius Rb/Rg 13 fm 6.5 fm 

average freeze-out proper time to 7.5 fm/c 

mean proper emission duration At 2 fm/c 

width of the space-time rapidity profile An 1.3 
pion mass 0.139 GeV/c^ 

kaon mass 0.493 GeV/c 2 



where the transverse expansion rapidity n t is assumed to scale linearly with the radial 
distance r: 

m = Vf-r— ■ (3i) 

' rms 

77/ is the scaling factor. Note that the transverse rms radius for the Gaussian source is 
r rms = V2R G , (32) 
while for the box-shaped transverse profile (29) one has 

r rms = ~-j= ■ (33) 

The same rms radius is thus ensured by setting Rb — 2Rq- The source dynamics shows 
up in the coupling of the velocity field (30) with the momentum (27) via the Boltzmann 
factor exp[— K ■ u/T}; T is the freeze-out temperature. Using (27) and (30) the scalar 
product in the exponent is written as 

K ■ u = M± cosh(F — i]) coshr/t — K± cos ip sinhr^ . (34) 

For our calculation we take the data-inspired parameter values listed in Table 1. 
They are not obtained from a careful fit to the data, but seem to reproduce the measured 
correlations from Pb+Pb collisions at 160 AGeV/c [27, 28] reasonably well. Note that 
the different values of r/f used in the two models lead to the same average transverse 
expansion velocity v± in both cases. This allows for a better comparison of the results, 
since v± appears to be the relevant (model-independent) physical quantity expressing 
the strength of the transverse collective expansion. The transverse sizes Rg and Rb 
are chosen such that the resulting rms radii in both models are equal. 

In order to provide the reader with a more intuitive picture we show in Fig. 1 
transverse cuts of these sources as seen by midrapidity pions of transverse momentum 
K± = 0.4 GcV/c. The different shape of the effective sources obtained from the two 
models is crucial for an understanding of the Mj_-dependences of the correlation radii. 
These will be studied in what follows. 
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Figure 1 Transverse cuts of the two models from Table 1 at 770 = for midrapidity pions 
(Y = 0) of transverse momentum K± = 0.4GeV/c. Left: Gaussian source; right: box-shaped 
source. (Note that the "corners" of the box-shaped source are a plotting artifact.) 



3.1 Af^-dependence of the correlation radii 

In Figs. 2 and 3 we show the Mj_-dependences of the correlation radii of all mentioned 
parametrizations. We evaluate the radii in the LCMS, and in order to get a non- 
vanishing cross-term and Yano-Koonin velocity we take particles at slightly forward 
rapidity (Y — 0, 770 = —1)- We compute both pion and kaon correlations. Instead of 
plotting the Yano-Koonin velocity in Fig. 3c we show the Yano-Koonin rapidity 

Y YK = lln\±^. (35) 

Z 1 YK 

In Fig. 3c' we show the modified YK rapidity, defined via an analogous relation involving 
v' . 

YK 

Let us comment on a few characteristic features. Although the average transverse ex- 
pansion velocity is the same in both models, the box-shaped source generates a stronger 
Mj_-dependence of R s than the Gaussian one. This effect is even stronger for R . More- 
over for higher transverse mass R is smaller than R s for the box-shaped source while 
the opposite is true for the Gaussian source. This results from the different shapes of the 
effective sources shown in Fig. 1: while for the Gaussian transverse profile the effective 
source for higher M± expands into the dilute tail of the density distribution, for the 
box-profile it is squeezed towards the outer boundary, leading to much smaller widths, 
especially in the outward direction. Note that this has important phenomenological 
consequences: even if the temperature is known, the exact amount of transverse flow 
cannot be inferred uniquely from the shape of R S (M±) without making assumptions 
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Figure 2 Correlation radii of the Cartesian parametrization for the models introduced in (26) 
with parameter values from Table 1. The calculation is performed in the LCMS for particle 
pairs at slightly forward rapidity (rjo = —1). Both pion (solid and dotted lines) and kaon 
(dashed and dash-doted lines) correlations were computed. Transverse profiles: box (solid, 
dashed), Gauss (dotted, dash-dotted). 



about the transverse source geometry. By combining measurements of R s and R 0l how- 
ever, one should be able to shed some light on this question. Definitely, this problem 
will deserve a more detailed study when fitting real data. 

The mentioned behaviour is also reflected in the M^-dependence of Rq. For a box- 
shaped source (x 2 ) < (y 2 ); via the last term in (20a) this leads to negative values of Rq. 
For kaons from this source the YKP parametrization does not even exist in the studied 
kinematic region. 

This last problem does not exist in the Modified YKP parametrization, as seen in 
Fig. 3. The behaviour of Rq is completely dominated by the last term in (24a). The 
comparison of both models is in accord with Fig. 2. 

Since the longitudinal geometry and dynamics are identical in both models, no 
differences are observed in the behaviour of the longitudinal correlation radii. Even in 




Figure 3 YKP (upper row) and Modified YKP (lower row) radii for the same models and 
the same kinematics as in Fig. 2. Note that R± = R' ± = R s - Line symbols are as in Fig. 2. At 
the given kinematics no YKP parametrization exists for kaons from the box-shaped source. 



case of i?'| 2 the correction terms (the last two terms in (24b)) lead to negligible effects, 
since they are multiplied by f3f which is small in the YK frame. 

3.2 YK rapidity and source velocity 

Since the direct access to the effective source velocity is probably the main strength 
of the YKP/Modified YKP parametrization, Fig. 4 is entirely devoted to the investi- 
gation of this particular feature. The velocity parameters are again encoded via the 
corresponding rapidities as introduced in (35). In Fig. 4a we show the characteris- 
tic dependence of the Yano-Koonin rapidity on the pair rapidity for a longitudinally 
(boost-invariantly) expanding source, as observed in the source center-of-mass (CMS) 
frame. The modified YK rapidity behaves in the same way. Note that the modified YK 
rapidity seems to be rather insensitive to the shape of the transverse density profile. 
This was already seen in Fig. 3c'. Fig. 4b (which is in fact a compilation of Figs. 3c 
and 3c') shows the Mj_-dependence of these rapidities for Y = 1 in the CMS. 

As discussed in the previous section, it is important to check to what extent the 
Yano-Koonin velocity really coincides with the velocity of the effective source. We also 
discussed the ambiguity of defining the latter for expanding effective sources. Here we 
study the difference Ydiff between the YK and/or Modified YK rapidity and the space- 
time rapidity (which is also the longitudinal fluid rapidity) of the point of maximal 
emissivity ij max (maximum of the emission function) . While for most studied cases this 
difference is positive for Y YK > and negative in the opposite case (i.e. the YK/Modified 
YK velocity is "faster" than the fluid velocity at the maximum emissivity point), it is 
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Figure 4 Investigation of the YK/YK' rapidities resulting from the two studied models: box- 
shaped (solid and dotted line) and Gaussian (dashed and dash-dotted) (see text). Quantities 
related to YKP parametrization are plotted by solid and dashed lines, those related to the 
YKP' one by dotted and dash-dotted lines, (a) YK/YK' rapidity as function of pair rapidity 
Y studied in the center-of-mass (CMS) frame of the fireball, K± = 0.3 GeV/c, note that dotted 
and dash-dotted line coincide; (b) the same quantities as functions of M±_ for Y = 1 calculated 
in the source CMS (r]o = 0); (c) the difference Yam between the YK/YK' rapidity and the 
space-time rapidity of the point of maximal emissivity (maximum of the emission function) 
as function of Y, K± — 0.3 GeV/c; (d) the same quantity but as function of M± for Y = 1, 
r)o = 0. 

the other way around for high M± pions from the box-shaped source. This difference is 
therefore again a model-dependent feature. However, it is important to stress that in all 
cases Ydiff is considerably smaller than Y" YK or Y^ K . Note also that, as discussed before, 
the rapidity ?7 max may well not be the "right" one to express the rapidity of effective 
source: if the source expands there is always an ambiguity of how to define its rapidity. 
Our choice is one of several posibilities. The small value of Ydiff should thus not be 
taken quantitatively as the difference between the experimentally accessible quantity 
Y" YK and the true source rapidity; its smallness rather indicates that the former maps 
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Figure 5 Correlation radii related to the emission duration (plotted by thick solid lines) as 
functions of M±. Upper row: box-shaped source, lower row: Gaussian source (see text). The 
remaining curves show the individual contributions: (a) and (b): /3^_{i 2 } (dashed), (x 2 — y 2 } 
(dotted) and — 2(3±{xt) (dash-dotted); (c) and (d): same as before but all terms divided by 
/3 2 ; (e) and (f): same as (c) and (d) but the dotted line corresponds to (x 2 )//3 2 . Units on the 
ordinate are fm 2 . 

the latter very well. We conclude that the measurement of longitudinal expansion in 
this way is reliable and does not depend on particular assumptions about the transverse 
geometry. This is valid for both YKP and YKP' parametrizations. 

3.3 The emission duration 

In order to analyze the measurement of the emission duration in more detail, we show in 
Fig. 5 the relevant correlation radii of all three parametrizations studied in this paper, 
together with the individual contributions out of which they are constructed. This 
study was performed for pions at Y = rjo = 0; in this case the LCMS and YK frame 
coincide and the decomposition is easy. This Figure can be summarized as follows: the 
geometric contribution (x 2 — y 2 ) is at least as important as the temporal contribution 
(t 2 ). This is even more true for the box-shaped source: here i? 2 < B? s and i? 2 < 0. The 
reason for this is the rapid decrease of the size in outward direction which was illustrated 
in Fig. 1. Note that such a behaviour was suggested in [25] as a possible indicator for 
opacity of the source, i.e., surface-dominated emission. We see here that already a 
source with a box-shaped transverse density profile exhibits this "opacity signature" for 
large K±, as clearly seen in Fig. 1. For us, however, it is important to conclude that, 
unless c 2 (t 2 ) is considerably larger than the geometric size of the source, the lifetime 
cannot be reliably extracted from interfcromctric measurement as one cannot really 
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correct for contributions from the transverse geometry. There exist many combinations 
of (t 2 ) and (x 2 — y 2 ) which all lead to the same interferometric signal. 

We finally note that, as shown in Figs. 5e,f, the radius parameter R'q in the modified 
YKP paramctrization does not allow for an extraction of the source lifetime, as already 
anticipated in Section 2. Rq is seen to be completely dominated by the term (x 2 )//3 2 
and practically insensitive to (t 2 ). 

4 Conclusions 

We have shown that the recently advocated YKP parametrization is not defined for 
all choices of kinematic variables. It is important to stress that this non-definiteness 
can only be identified experimentally through comparison of fit results from different 
parametrizations via the relations introduced in Appendix B. As it is hard to interpret 
fit results from a possibly ill-defined paramctrization it is absolutely necessary to per- 
form such comparisons and cross-checks, and one should worry seriously about cases 
where inconsistencies appear! 

A Modified YKP parametrization was introduced which does not suffer from these 
problems. Eventually, it therefore also can serve for the check of the YKP radius 
parameters extracted from a fit to the data [27]. The corresponding relations were 
given in Appendix B. 

It was shown that the velocity parameters of both YKP and Modified YKP parame- 
trizations accurately reflect the longitudinal velocity of the effective source. This feature 
appears to be fairly model- independent. The idea [14] of introducing these velocity 
parameters as the effective velocity of a non-expanding effective source thus remains 
valid in the context of strongly expanding emitters. 

Our consideration of the measurement of the emission duration showed that it is 
practically impossible to extract this quantity in a model-independent way. On the 
other hand, clear differences in the correlation patterns resulting from models with 
Gaussian and with box-shaped transverse profiles were observed. These could help to 
decide between models which are more or less suitable for a description of the data. 
In fact, recent studies show that a box-shaped source fits the data of the NA49 collab- 
oration better than a Gaussian one [28]. By making such specific assumptions about 
the transverse source geometry (and only then) the emission duration can be extracted 
from the data fit. 
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A Definition range for the YKP/Modified YKP parameters 

The Yano-Koonin velocity v YK is only defined if the discriminant 

D = (A + B) 2 — 4C 2 , (36) 

is positive. Here we study the conditions for that. From (24) it is clear that 

A' > , (37a) 
B' > , (37b) 

so we conclude 

\A' + B'\ = A' + B' . (38) 
It can be proven that 

D' = {A' + B'f - 4C" 2 > , (39) 

or 

\A' + B'\ -2\C'\ > 0, (40) 

i.e., the modified Yano-Koonin velocity is defined everywhere. Indeed, due to (38) 
inequality (40) may be written as 

A' + B't2C">0, (41) 

where the upper sign stands for C > and the lower for C < 0. Inserting expressions 
(24) leads to 

^'^({(i-^M.-Ax)}^. (42) 

This proves (39) - the modified YKP parametrization is defined everywhere (except, of 
course, the point K± = j3±_ = [12, 16]). 

From expressions (24) one also can see that 

A + B = A' + B'-^-(f), (43) 



1C = 2^-||(y 2 >. (44) 



It is therefore possible to obtain negative values for D: 

D=(A + B) 2 - 4C 2 < . (45) 
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For example, if fli — and (y 2 ) has a value bigger than (x 2 ), there exist values for f3±_ 
such that 

\A + B\ = \A> + B' — ^(f}\ < 2\C'\ = 2\C\ . (46) 

Hence, there are kinematic regions where the YKP parametrization is not defined. (For 
example see Section 3.) 



B Relations between correlation radii of different parametrizations 

In this Appendix relations between correlation radii of various parametrizations are 
listed. 

The simplest and most obvious one is 

R 2 s =Rl=R'l. (47) 
The remaining YKP parameters are related to the modified ones in a simple way: 

(jfl+t^)- * jtf, (48a) 



y 2 



B = ^(jtf + ^jtf)-^, (48b) 

C = l' 2 < K (K 2 + R?)-^R'L (48c) 
and the inverse relations are given by 

A' = 7 2 (i? 2 + W 2 KJ R 2 ) (49a) 

B' = 7 2 (i? 2 +«)+|M?i, (49b) 

C = 7 2 VYK (i? 2 + i? 2 )+|-i?i. (49c) 

Similar relations exist, of course, between the original and/or modified YKP param- 
eters and the Cartesian correlation radii. The latter are expressed in terms of YKP 
parametrization [12, 16] 

Rl-R 2 s ^ R\ m = (r\ + v 2 YK Rj ) , (50a) 

Rf = {l-^Rl+^ifr-v^Y^Rl + Rf) , (50b) 

Rli = /3x(-Ai? 2 +7 2 (A-« YK )(^o + ^ 2 )) , (50c) 
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while the inverse relations read 
1 

Pi 

Pi R 2 , Pf 



A = ^R 2 diS , (51a) 
B = Rf-2^Rl l + ^R 2 diS , (51b) 



C = -^R 2 oi + ^R 2 diS - (51c) 

If instead of the original one the Modified YKP parametrization is used then 

Rl = Ph' 2 (R'o 2 +v' Y 2 K Rf), (52a) 

Rf = (1-A 2 )i?f+7' 2 (A-< K ) 2 (<+^| 2 ) , (52b) 

Rl = Pi- (-PiR\ 2 + 7 ,2 (A - <J (< + if)) , (52c) 



and 



1 



A' = j^Rl, (53a) 
B' = Rf-2^-Rl l + ^-Rl, (53b) 

c = -^Rli + w R2 o- ( 53c ) 



Note that the only formal change in (52)/(53) if compared with (50)/(51) is the re- 
placement of i?^ iff with Rl and the use of primed radius parameters. 

C Off-shell corrections to the YKP/Modified YKP correlation radii 

In this Appendix we list the correction terms to the Gaussian width parameters of the 
correlator which arise from releasing the approximations made in (4), namely: 

(Ad) setting Pi{pi)Pi{p2) — * [Pi{K)\ 2 m denominator of that relation - these correc- 
tions are denoted with 8d\ 

(An) replacing K with its on-shell value in the emission function appearing in the 
numerator, i.e., K° — ► E K = \J K 2 + m 2 - corresponding terms will be assigned 
with S n . 

Here 

Pi(p) = f d^xSix^p), (54) 

denotes the single-particle invariant momentum distribution. In fact, only the correction 
terms to the shorthands (20) and (24) are derived; those for the YKP/Modified YKP 
radii can then be obtained via (19)/(23). 
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The following relations are obtained by combining the correction terms known for 
the Cartesian parametrization [17]. They read for (Ad) 



5 d Rl 



S d R 2 = 
S d Rf = 
SdRli = 

and for (An) 

S n R 2 = 
S n Rl = 



S n Rf 



4K ± dK ± 
1 d 2 



In Pi (if) 



In Pi (if). 



4dK 2 ± 



1 



d 2 



4 dK ± dKi 



In Pi (if), 



1 d 



2 dm 2 



In Pi (if) , 



5nR\ 



ol 



(55a) 
(55b) 
(55c) 
(55d) 

(56a) 
(56b) 
(56c) 
(56d) 



Now we have to combine these terms according to (51) and (53). The most trivial 
correction, of course, is the one belonging to R\ and R' 2 as these radius parameters 
arc identical with R 2 : 



5dR±_ — 5dR'±_ — SdR s — 



1 



4if_L dK ± 
For the remaining parameters one finds 



In P (if). 



S d A 
S d B 
S d C 



d 2 



dK\ K± dK\ 



V 2 



0? 1 



1 _ d 

o 
d 



V + 



pj_ K ± dK ± 
0± K ± 



dK ± 



In Pi (if). 



In Pi (if) . 



where two operators have formally been introduced 



V = 
V 2 = 



Pi d 



d 

dKi (3±_ dK ± ' 
d 2 2ft d 2 



Pf d 2 



dKf p± dKidK ± Pi dKi 



2 ' 



(57) 

(58a) 
(58b) 
(58c) 

(59a) 
(59b) 
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Corresponding relations for the modified radii read 



S d A' 
S d B' 

s d c 



1 



\v 2 lnfMtf), 



1 1 ^ d 
V 

A 13 ± dtf± 



InPi(A-). 



(60a) 
(60b) 
(60c) 



On the first sight, these correction terms diverge as (3± 



gence cannot be a priori excluded in (20)/(24). One might argue that (x 2 — y 2 



0. In fact, such a diver- 

(tx), 

and (zx) in (20) should vanish for (i±_ — > due to restored azimuthal symmetry; it is, 
however, not guaranteed that they vanish sufficiently fast. Fortunately, the apparent 
singularities actually cancel: rewriting M± = Kj_ + m 2 



dK ± 
d 2 



K ± d 

Ml dM ± ' 

m 2 d K 2 
+ 



M\ dM ± 



we find 
8 d A 



El 



AMI 



Ml dM\ ' 



1 



dM 2 



S d B = 



M± dM ± 
d 



InPi(A-), 



5 d C = - 



1 E K 

A Mi 



dM ± 



V 



Ml 



InPi(A-) 
d 



In Pi (iT) 



and 



4M 2 



d 2 1 to 2 d 

dM 2 + 0[ E 2 K M± dM±_ 



InP^K). 



—-2 

X" + 



1 E K 

~4 Mj_ 



Pi Ml dM ± 



~ ~p\ E K Ml dM ± 



In Pi (if), 
d 



lnPi(K). 



with 



D' = 



V 



d 

dKi 
d 2 
dK 2 



Pi E K d 



Mi dM ± ' 
2B,E K d 2 



M± dKi dM± 



+ 



PfE 



K 



M 



dM 2 



(61a) 
(61b) 

(62a) 
(62b) 
(62c) 

(63a) 
(63b) 
(63c) 

(64a) 
(64b) 



From (62) one sees that the correction terms belonging to the YKP parametrization 
behave well. The singularities cancel thanks to the appearance of P^ iff in (51). In 
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case of the modified parametrization R 2 appears instead of -R^iff m (53), and thus the 
correction terms (63) still diverge for vanishing f3±. 

Corrections due to approximation (An) are obtained in a similar way as above. The 
one corresponding to R\ (and R' 2 as well) is again obvious 

' "' lnPi(if), 



<5nPj_ — S n R' 2 — S n R 2 s 
KP i 
In Pi (if), 



2 dm 2 

and the relations for YKP parametrization are 
1 d 



S n A 

S n B 
S n C 



2 dm 2 
1 d 



2 dm 2 



In Pi (if), 



0. 



For the modified YKP parameters one finds 
1 d 



S n A' 2 
S„B' = - 

S n C = - 



1 + 



dm 2 



In Pi (if), 



1 

2 

1 0i d 

2 2 dm 2 



0\) dm 2 



In Pi (if), 



In Pi (if). 



(65) 



(66a) 

(66b) 
(66c) 



(67a) 
(67b) 
(67c) 



While the corrections to the YKP parameters are clearly regular, this is not true for 
the modified ones. These may be rewritten by using 

d 1 J d_ 

dm 2 ~ 2 M± dM ± ' ( ' 

This leads to 

1 1 0i-l 



SnA' 
SnB' 



4 M_l 0\ dM ± 



In Pi (if), 



1 



A 2 



i 



- 1 + -t) lnPi(if) , 

4 V 0\) M ± dM ± u 

1 



*» c ' = -mwidk M{K) - 



(69a) 
(69b) 
(69c) 



If now these corrections are added to those due to (Ad) displayed in (63), one finally 
arrives at the total correction terms to the modified YKP radii: 



SA' 
SB' 
SC 



1 1 

4 Ml 



E\ d 2 



V - 

1 E K 

' 4 M± 



M± dM\ 
1 



f 



F? 
K 



d 



dM\ 



In Pi (if). 



V 



1 

Pi E K 
Ml 



013. 

Ml 



dM ± 



In Pi (if), 



dM ± 



In Pi (if). 



(70a) 
(70b) 
(70c) 
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Note that no divergences survive in the final expressions. 
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